Statistical correlations in an ideal gas of particles obeying fractional exclusion 
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After a brief discussion of the concepts of fractional exchange and fractional exclusion statistics, 
we report partly analytical and partly numerical results on thermodynamic properties of assemblies 
of particles obeying fractional exclusion statistics. The effect of dimensionality is one focal point, 
the ratio n/k&T of chemical potential to thermal energy being obtained numerically as a function 
of a scaled particle density. Pair correlation functions are also presented as a function of the 
statistical parameter, with 'Friedel' oscillations developing close to the fermion limit, for sufficiently 
large density. 
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I. INTRODUCTION 

Current interest in fractional statistics pj, 0, HJ is mo- 
tivated by its possible relevance for fractional quantum 
Hall effect 0, 0| and high-temperature superconductiv- 
ity ((J. In this context, an early application was made by 
Lea et al. 0, @| to discuss the de Haas-van Alphen oscil- 
latory orbital magnetism of a two-dimensional electron 
gas (2DEG). Quite recently, it has been proposed that 
noise experiments in quantum Hall fluids might reveal 
the existence of elementary excitations obeying fractional 
statistics (9l. Hoj|. 

Fractional exchange statistics [l|, [H, H[ arises when the 
many-body wave function of a system of indistinguish- 
able particles (dubbed anyons) is allowed to acquire an 
arbitrary phase e M7r upon an adiabatic exchange pro- 
cess of two particles. Here, a is the so-called statistical 
parameter, interpolating between a — (bosons) and 
a = 1 (fermions). Such an exchange produces a non- 
trivial phase only if the configuration space of the col- 
lection of particles under study possesses a multiply con- 
nected topological structure. Therefore, fractional ex- 
change statistics is usually restricted to two spatial di- 
mensions, d = 2 (see Ref. [ll[ for a review). However, 
fractional exchange statistics can be formalized, to some 
extent, also in d = 1 (l2l. Il3j. 

A different concept of fractional statistics, namely frac- 
tional exclusion statistics, is based on the structure of 
the Hilbert space, rather than the configuration space, of 
the particle assembly, and is therefore not restricted to 



d < 2. Fractional exclusion statistics has been introduced 
by Haldane [14j . who considered the ratio 
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-AD/AN, 



(1) 
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where AD denotes the change in size of the subset of 
available states in the Hilbert space corresponding to a 
change A A of the number of particles (i.e. elementary 
excitations). Clearly, one has again g — for bosons 
and g = 1 for fermions, the latter being a consequence of 
Pauli exclusion principle. At variance with anyons, par- 
ticles obeying fractional exclusion statistics are usually 
dubbed g-ons or exclusons. The distribution function for 
fractional exclusion statistics has been derived by Wu 
[l5| . A limiting form of the same distribution function 
had been derived by March et al. [1 €31 ] within an a pprq x- 
imate chemical collision model (see also Refs. [171 . Ha|). 
Following Wu's distribution function [l5[ , the thermody- 
namics of an ideal gas with fractional exclusion statistics 
has been studied in some detail in arbitrary dimensions 

[SH. 

The relation between fractional exchange and exclu- 
sion statistics is elusive. In particular, it has been em- 
phasized [2l[ that in order to derive a consistent statisti- 
cal mechanics for g-ons, Haldane's generalized exclusion 
principle, Eq. |T]), must hold locally in phase space, which 
indeed applies rigorously to either bosons or fermions. In 
other words, AD in Eq. |T]) should be related to states of 
close energy, which is brought about by an effective in- 
teraction which is local in momentum space. From this, 
it has been concluded that anyons are not ideal g-ons, 
but interacting g-ons [2l| . Haldane 14 1 provides explicit 
examples of systems for which fractional exclusion, albeit 
not exchange, statistics might be applicable. An earlier 
attempt to relate the parameter g of fractional exclusion 
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statistics and the statistical parameter a of fractional ex- 
change statistics has been made in Ref . [22[ , on the basis 
of the virial expansion of the g-ons partition function. 
More recently, an analytic, monotonic relation g = g(a) 
has been derived in Ref. [23j in the case g = 1/m (for 
integer m), which is relevant for the fractional quantum 
Hall effect (see also Ref. [24j for a review). 

In an electron gas, statistical correlations between non- 
interacting fermions arise as a consequence of particle 
indistinguishability. In the case of fermions, the anti- 
symmetry of the wave function implies that two such 
particles can never be found simultaneously at the same 
point in space. As a result, each electron is surrounded 
by an 'exchange hole', a region in which the density of 
same-spin electrons is smaller than average and in which, 
therefore, the positive background charge is not exactly 
cancelled. It is the interaction of each electron with the 
positive charge of the exchange hole that gives rise to the 
exchange energy. The situation is quite the opposite for 
bosons, where statistical correlations are attractive and 
may be thought of being ultimately responsible of Bose- 
Einstein condensation. In the case of intermediate or 
fractional statistics, the question naturally arises of how 
many particles does a particle effectively 'see' around it- 
self. In other words, one is led to consider the effect of 
Haldane's generalized exclusion principle on the exchange 
hole which a <?-on 'digs' around itself. 

In this work, we attempt to answer such a question, 
by evaluating a suitably defined generalization of the pair 
correlation function between g-ons (i.e. particles obeying 
fractional exclusion statistics), in arbitrary dimensions d 
(d = 1, 2, 3), as a function of temperature T and particle 
density. Earlier results along this direction include the 
seminal work of Sutherland [25| , for a system of fermions 
or bosons in d = 1 interacting through a singular poten- 
tial ~ 1/f 2 , a problem which can be mapped into anyons 
[ll| . The evaluation of a pairwise correlation function for 
anyons (i.e. particles obeying fractional exchange statis- 
tics) in d = 2 at T = has been also recently analyzed 
in Ref. 26], in connection with an intensity interferome- 
try gedanken experiment, which may prove useful in the 
context of quantum computing. The question is also in- 
timately related to that of finding Uhlenbeck's statistical 
interparticle potential [27[, which has been studied for 
the two-anyons system in d = 2 [28| • 

The paper is organized as follows. In Sec. [IT] we derive 
the pair correlation function for a homogeneous assem- 
bly of g-ons as a function of interparticle separation r. 
In Sec. Ill Al we study numerically the dependence of the 
chemical potential on density at a given temperature, 
for arbitrary dimensions d = 1,2,3, thereby reproducing 
Wu's analytical result at d = 2. In Sec. Ill Bl we study nu- 
merically the pair correlation function, and in particular 
focus on the possible occurrence of 'Friedel' oscillations, 
depending on the statistical parameter and on the scaled 
density. We eventually summarize in Sec. Mil 



II. PAIR CORRELATION FUNCTION FOR 
NON-INTERACTING p-ONS 

Statistical correlations between indistinguishable, 
hard-core particles can be described by means of the pair 
correlation function g(r), defined as the normalized prob- 
ability of simultaneously finding a particle at position r 
and a particle at position r = 0. (Here and below, we 
are implicitly neglecting spin.) In the formalism of sec- 
ond quantization, the pair correlation function can be 
written as 



9(r) = 



(#t( r )$t(o)#(o)#(r)) 
n(r)n(0) 



(2) 



Here, vF^r) [^(r)] is a creation [annihilation] quantum 
field operator at position r, n(r) = (^(r)^^)) is the 
single-particle probability density at position r, and (. . .) 
denotes a quantum statistical average associated with the 
equilibrium distribution of the particle assembly under 
study (see below) . In the case of a homogeneous system, 
one obviously has n(r) = n(0) and g(r) = g(r). 

In the case of either fermions or bosons, the four-point 
average appearing in the right-hand side of Eq. ([2]) can 
be reduced by means of Wick's theorem as 



(*t( r )^,t( )^(0)*(r)) 



(* t (r)*(r))(* t (0)*(0)) 
±(* t (r)*(0))(* t (0)*(r)), 

(3) 



the plus (minus) sign corresponding to bosons (fermions). 
In the case of particles obeying exclusion statistics, the 
appropriate generalization of Wick's theorem depends on 
the behavior of excluson field operators under exchange. 
As emphasized in the Introduction, the relation between 
fractional exclusion and exchange statistics is not com- 
pletely settled. Within fractional exchange statistics in 
reduced dimensionality (d < 2), one obviously has the 
following graded commutation relations, 

[*(r),*(r')]a = *(r)*(r') -e M7rs ( r ' r ')*(r')*(r) 

= 0, (4a) 

[*t(r),*t(r% = 0, (4b) 
[*(r),*t (r % = S(r-r>), (4c) 

with s(r,r') = ±1, depending on whether the exchange 
is being performed along a (counter)clockwise path in 
d = 2 [nj], or s(r,r') = sgn(r-r') ind = 1 [H. Eqs. Q 
naturally reduce to the familiar commutation and anti- 
commutation relations for bosons (a = 0) and fermions 
(a = 1), respectively. In d = 3, the sign s(r,r') is ill- 
defined, and in fact fractional exchange statistics does 
not apply to dimensions d > 3. In order to have gener- 
alized commutation relations valid for exclusons in arbi- 
trary dimensions, we therefore assume real field opera- 
tors, assign a statistical exchange parameter a to exclu- 
sons (g-ons), and postulate that standard quantum Bose 
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or Fermi statistics are weakly violated as in 



[*(r),*(r% 



*(r)*(r') - cos(a7r)*(r')^ r (r) 
0, (5a) 
0, (5b) 
S(r-r'), (5c) 



which again reduce to the familiar relations for bosons 
and fermions, in the appropriate limits. Following 
Ref. (29[ for the generalization of Wick's theorem cor- 
responding to Eqs. ©, Eq. ([3]) becomes 

(*t( r )^t( )*(0)W(r)) = (* t (r)*(r))(* t (0)*(0)) 
+ cos(mr) (* t (r)*(0)) (* f (0)*(r)). 

(6) 

Correspondingly, the pair correlation function for a ho- 
mogeneous assembly of g-ons reads 



g(r) = 1 + cos(a7r) 



|(ttt(r)*E(0))| s 
n 2 {0) 



(7) 



which correctly reduces to the fermion limit (a = 1), with 
g(r) < 1, manifestly [30l |. 

For translationally invariant systems, one may use a 
plane wave expansion to find 



/(r) = 1 + cos(a7r) 



n(r) 



n(0) 



where 



n(r) 



d d k 
(2^p 



r n(k) 



(8) 



(9) 



is the Fourier transform in d dimensions of the single- 
particle distribution function n(k) for g-ons in equilib- 
rium at temperature T. 

Following Wu [l5j |. the average occupation number n, 
of exclusons in an energy state Ci at equilibrium is 



1 



w(0 



(10) 



where £ = exp[(e^ — /z)/A;bT], /x is the chemical potential, 
and the function w(Q satisfies the functional equation 



w a (l 



(11) 



Here and below, we shall use of the standard notation 
g = a, even for exclusons. Eq. (fT0|) establishes the ther- 
modynamics of the excluson assembly, once the energy 
level distribution tj is given, which is determined by the 
free-particle dynamics. In order to make contact with 
the ordinary Bose and Fermi limits, we shall then as- 
sume that €i = h 2 k 2 /(2m), and then treat kj i— > k as 
a continuous variable, as is understood in Eq. (jSJ) with 
7ii i— ► n(k). 

An approximate form of the distribution function for 
exclusons had been earlier derived by March et al. [HI, 



[ijiLLll from collision theory, by using the detailed balance 
hypothesis. Their result was that n^ 1 — £ + a, with a 
interpolating between Bose (a = —1) and Fermi (a = +1) 
statistics. Comparison with Wu's result, Eq. (|10p. then 
yields [U, Qj| ' 



a(a, w) = w + a — w a (l + w) 1 



2a -1, (12) 



the latter approximation, independent of w, holding in 
the limit w ^> 1. Eq. (fT2)) already correctly interpolates 
between the Bose (a = 0, a = — 1) and Fermi (a = 
1, a = 1) limits. Of course, as to becomes small, the 
dependence of a on w is significant, though 2a — 1 is not 
a bad approximation even at w = 1 (giving for a 
a = instead of 0.09, and for a — h, a 
-0.44). 



2 ' 

i rather than 



A. Chemical potential 

Within these assumptions, the general properties and 
some exact results of the thermodynamics of an assembly 
of noninteracting exclusons have been derived to some ex- 
tent [15|, [l9j, l20( . Some properties of Wu's function w(C) 
in the complex £-plane and their relevance for Cooper 
pairing of exclusons have been discussed in Ref. [3l|. In 
particular, by evaluating the particle density 



N 
V 



J2^y 



n(k) 



(13) 



in dimensions d = 2, Wu [l5j was able to derive the re- 
lation between chemical potential fi, temperature T and 
particle density N/V explicitly as 



X 2 N 

i 

V 



log 



1 — exp 



\ 2 N 



(14) 



where A = hyj 2n / {mk^,T) is the thermal wavelength. 
Performing the integration in Eq. (|13| for a homogenous 
system in arbitrary dimensions d (see also Ref. [20j), such 
a relation can be generalized implicitly as 



dw — j- r( 

wo w(l+w) 



(15) 



where T(x) is Euler's function, wq = wo(fj,/ k^T) is a 
generalized inverse fugacity, implicitly defined by Wu's 
functional equation for e = 0, i.e. 



and 



log 



<(l + u>o) 



w a (l + w) 1 - a e^ kBT 



,-H/kB.T 



(16) 



log 



w a {l 



w 



w^l + wo) 1 " 



Eliminating wq between Eqs. (|15p and (|16p. we obtain 
the relation between the ratio /x/fceT of chemical poten- 
tial to thermal energy and the scaled density X d N/V, 
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X d N/V 



FIG. 1: Chemical potential vs scaled density, Eq. (|15p . for 
a = 0, 0.5, 1 (bottom to top) and d = 1, 2, 3. 

parametrically. This is plotted in Fig. [1] for several val- 
ues of the statistical parameter a, and d = 1,2,3. We 
explicitly note that the dilute limit \ d N/V <C 1 corre- 
sponds to w > 1 (fi/ksT — > — oo). We thus recover a 
monotonic relation, as expected, which for a — (Bose 
limit) leads to an infmitesimally small negative chemical 
potential, in the limit of large density. 

B. Pair correlation function 

Following an analogous procedure, and performing the 
integration required in Eq. ^ for an isotropic system in d 
dimensions, for the Fourier transform of the equilibrium 
distribution function of an assembly of noninteracting ex- 
clusons we find 

X d n(x) = x~ u (°° dw — — ^ Tt u/2 J V hxl 1 ' 2 ) , (18) 

J Wo w(l+w) V J 

where x = y^r/A, v = (d— 2)/2, and J v (z) is a Bcsscl 
function of first kind, arising from the integration over 
angles in d dimensions. Making use of the asymptotic 
properties of the Bessel functions [Hj], it can be shown 
that Eq. ([TH]) correctly reduces to Eq. (T5")) for the particle 
density in the limit x — > 0. Further analytical results 
concerning the local excluson density n(x) are derived in 

App.m 

Figure [3] shows g(x) in d = 2 for u; = 0.1 (corre- 
sponding to a relatively large particle density), and for 
a = — 1. In the fermion limit (a — 1) and close to it 
(0.5 < a < 1), the pair correlation function exhibits a 
correlation 'hole' around x = 0, whose depth decreases 
with decreasing a and eventually vanishes as a — > \. 
In the same range of values of the statistical parameter, 
g(x) is characterized by damped 'Friedel' oscillations (cf. 
inset of Fig. which are more pronounced close to the 
fermion limit. These oscillations are absent within nu- 
merically accuracy in the prevalently bosonic range of 
the statistical parameter, < a < 0.5, where the pair 
correlation function displays a monotonically decreasing 
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x 



FIG. 2: Pair correlation function g(x) (in scaled units), 
Eq. JSJ, in d = 2, for wo =0.1 (dense limit) and a = — 1 
(top to bottom). Inset shows Friedel oscillations, which are 
present only close to the FD limit (0.5 < a < 1). 
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0.5 1 1.5 2 2.5 



x 

FIG. 3: Pair correlation function g(x) (in scaled units), 
Eq. JS), in d = 2, for Wo = 0.5 (dilute limit) and a — 0—1 (top 
to bottom). Inset shows the absence of Friedel oscillations, 
even close to the FD limit (0.5 < a < 1). 



behavior, with correlations 'piling up' at x = 0, in con- 
trast with the fermion limit. Qualitatively similar results 
(not shown here) are obtained in d — 1 and d = 3. 

Figure [3] shows again the pair correlation function in 
d = 2 for a = — 1, but now for wq — 0.5, correspond- 
ing to a lower particle density. While the fermion side 
(0.5 < a < 1) is still characterized by a correlation hole, 
and the boson side displays a monotonic decrease of cor- 
relations, with a maximum at x — 0, 'Friedel' oscillations 
are now beyond graphical resolution even in the quasi- 
fermion regime (cf. inset of Fig. [3]), as an effect of the 
reduced particle density. Again, qualitatively similar re- 
sults obtain also in d — 1 and d = 3. 

In order to make more quantitative the latter state- 
ment, we now inquire on the condition giving rise to 
'Friedel' oscillations in the pair correlation function. 
Generally speaking, Friedel oscillations in many-body 
properties of fermion assemblies arise from the presence 
of a discontinuity in the Fermi distribution function at 
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FIG. 4: Phase diagram in the plane of scaled density X d N/V 
vs statistical parameter a. The phase diagram features re- 
gions where pronounced 'Friedel' oscillations occur in the pair 
correlation function (F), and where such oscillations are ab- 
sent, or beyond graphical resolution (nF). Different lines refer 
to various dimensionalities (d — 1, 2, 3). 

the Fermi level. Besides decreasing at T = as a function 
of distance with a power law, depending on dimensional- 
ity, such oscillations are then exponentially damped with 
a characteristic decay length ~ A 2 fcp, with fcp the Fermi 
wavevector. 

In the case of particles obeying fractional exclusion 
statistics, Wu's distribution function Eq. (jTOJ) in scaled 
units displays a smoother inflection point at a value e* 
corresponding to w* — a — 1 + \/a 2 — a + 1 in Eq. (fl~7f . 
When u>» > wq, such an inflection point lies within the 
integration range in Eq. ([18]), and does therefore give rise 
to marked 'Friedel' oscillations in the behavior of the pair 
correlation function. Such a condition has been stud- 
ied numerically as a function of scaled particle density 
\ d N/V and statistical parameter a, and is reproduced 
in Fig. 2] for various dimensionalities d = 1,2,3 in the 
form of a phase diagram. 

The plane of scaled density X d N/V vs statistical pa- 
rameter a is divided into two regions by the line implic- 
itly defined by the equation w* = wq. At a fixed value 
of the statistical parameter a, say, one might think of 
increasing the fermionic-like correlations of a given g-on 
assembly (i.e. sharpening the 'Friedel' oscillations in the 
pair correlation function) by suitably increasing its scaled 
density X d N/V, i.e. increasing density N/V or reducing 
temperature T. 

Within this picture, one may also estimate the wave- 
length A of the oscillations in the pair correlation func- 
tion with the approximate condition 2y / e7A « 2ir, which 
straightforwardly leads to 



where all dimensions have been restored. Clearly, as 
e* — > as a — > 0, one recovers the absence of 'Friedel' 
oscillations in the bosonic limit. 



III. SUMMARY 

A brief discussion of the different concepts of fractional 
exchange and fractional exclusion statistics is first given. 
Then we have focussed on some thermodynamical prop- 
erties and on the pair correlation function of an assembly 
of non-interacting particles obeying exclusion statistics in 
arbitrary dimensionality d = 1,2,3 and given tempera- 
ture T. In particular, Eqs. (fT5|) and (fl"6|) allow the chemi- 
cal potential to be studied as a function of scaled particle 
density and dimensionality. The pair correlation function 
is then studied as a function of interparticle distance for 
different values of the statistical parameter a and of the 
scaled density X d N/V, where A is the de Broglie thermal 
wavelength. The behavior of the pair correlation func- 
tion correctly exhibits a correlation 'hole' in the fcrmion 
limit, and a pronounced maximum in the boson limit. 
Away from the boson limit, it also features 'Friedel' os- 
cillations, whose sharpness and wavelength can be tuned 
as a function of the statistical parameter and of scaled 
particle density. Our results show that, in principle, tun- 
ing density and temperature may sharpen the statistical 
correlations among non-interacting exclusons, thus end- 
ing in a behavior closer to cither the fermionic or the 
bosonic limit, depending on dimensionality and on the 
value of the statistical parameter a. This may be rele- 
vant to systems exhibiting the fractional quantum Hall 
effect flij . where anyons, i.e. particles obeying fractional 
exchange statistics, can be regarded as effectively inter- 
acting exclusons [2l| . 
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APPENDIX A: RECURSION FORMULAS FOR 
THE LOCAL EXCLUSON DENSITY 

Eq. (|18[) yields n v (x, wq) = n(x) as a function of scaled 
interparticle distance x and generalized inverse fugacity 
wq, at a given statistical parameter a and reduced dimen- 
sionality v = (d — 2)/2. By differentiating Eq. (fl~5|) with 
respect to x, and making use of the appropriate recursion 
formula for the Bessel functions [32| , one obtains 

^ + 2xh„ +1 = 0. (Al) 

Similarly, by differentiating the same Eq. (JTSJ) with re- 
spect to wq, and taking into account Eq. (jTTJ) , one ob- 
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tains 



where 



and 



dw 



dw 



n v -i + p v (w ) = 0, 



a + w 



dwj w (l + w )' 



(A2) 



(A3) 



cursion formulas relating the local excluson density at dif- 
ferent dimensionalities. Combining Eqs. (|Alj) and 1|A2|) . 



and noting that p„ + i(w;o) 
d = 1,2, 3, one obtains 



lim„ 



t p v+ i(w) = for 



d 2 h v 
dwodx 



= 2x 



W 



wq(1 + w ) 



(A5) 



Pu{w) 



r(i/+ 1) w(i + w) 



(A4) 



is the density of states in the w variable, in d = 2v + 2 
dimensions. Eqs. (|Alj) and (|A2|) are then differential re- 



which is a hyperbolic partial differential equation for the 
local excluson density, similar to Klein-Gordon equation, 
but with a variable 'mass' term. 
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